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ABSTRACT 

We measure the radio luminosity function (RLF) of steep-spectrum radio sources us- 
ing three redshift surveys of flux-limited samples selected at low (151 & 178 MHz) 
radio frequency, low-frequency source counts and the local RLF. The redshift surveys 
used are the new 7C Redshift Survey (7CRS) and the brighter 3CRR and 6CE surveys 
totalling 356 sources with virtually complete redshift z information. This yields un- 
precedented coverage of the radio luminosity versus z plane for steep-spectrum sources, 
and hence the most accurate measurements of the steep-spectrum RLF yet made. We 
find that a simple dual-population model for the RLF fits the data well, requiring 
differential density evolution (with z) for the two populations. The low-luminosity 
population can be associated with radio galaxies with weak emission lines, and in- 
cludes sources with both FRI and FRII radio structures; its comoving space density p 
rises by about one dex between z ~ and z <~ 1 but cannot yet be meaningfully con- 
strained at higher redshifts. The high-luminosity population can be associated with 
radio galaxies and quasars with strong emission lines, and consists almost exclusively 
of sources with FRII radio structure; its p rises by nearly three dex between z ~ and 
z ~ 2. These results mirror the situation seen in X-ray and optically-selected samples 
of AGN where: (i) low luminosity objects exhibit a gradual rise in p with z which 
crudely matches the rises seen in the rates of global star formation and galaxy merg- 
ers; and (ii) the density of high luminosity objects rises much more dramatically. The 
integrated radio luminosity density of the combination of the two populations is con- 
trolled by the value of p at the low-luminosity end of the RLF of the high-luminosity 
population, a quantity which has been directly measured at z ~ 1 by the 7CRS. We 
argue that robust determination of this quantity at higher redshifts requires a new 
redshift survey based on a large (~ 1000 source) sample about five times fainter than 
the 7CRS. 

Key words: radio continuum: galaxies - galaxies: active - quasars: general - 
galaxies: evolution 



1 INTRODUCTION 

The radio luminosity function (RLF) seeks to derive from 
observed samples and surveys of radio sources, their space 
density per unit co-moving volume and how this changes 
with source luminosity. It is an essential input to both grav- 
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itational lensing studies which probe cosmic geometry (e.g. 
Kochanek 1996) and modelling of the clustering of radio 
sources (Magliocchetti et al. 1998) . The RLF derived at low 
radio frequencies is important for understanding the con- 
tent of high frequency surveys and hence deriving jet beam- 
ing parameters (e.g. Jackson & Wall 1999). The shape and 
evolution of the RLF provide important constraints on the 
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nature of radio activity in massive galaxies and its cosmic 
evolution. 

Longair (1966) was one of the first to attempt to deter- 
mine the evolution of the radio source population. At that 
time there were very few radio source redshifts known and 
the main constraint upon his models came from the low- 
frequency source counts down to S151 w 0.25 Jy. Longair 
found that the data were best-fit by models where the most 
powerful radio sources undergo greater cosmic evolution (in 
comoving space density p) than less-powerful sources. The 
subsequent acquisition of a substantial fraction of the red- 
shifts for the 3CRR complete sample of Laing, Riley & Lon- 
gair (1983) enabled the evolution to be better constrained. 
Using the U/ V max test (see Section ^) they concluded that 
the most powerful radio sources undergo substantial evolu- 
tion with similar evolution for both powerful radio galaxies 
and quasars. Unfortunately, the tight correlation between ra- 
dio luminosity and redshift in this sample meant that they 
were unable to make any progress on the form of the evolv- 
ing RLF. Wall, Pearson & Longair (1980) suggested that 
complete samples a factor of ~ 50 times fainter than the 
3CRR sample would be necessary to differentiate between 
possible models of the RLF. 

The most comprehensive study of the radio luminos- 
ity function prior to this work was performed by Dunlop 
& Peacock (1990; hereafter DP90). Using several complete 
samples selected at 2.7 GHz with lower flux-limits rang- 
ing from 2.0 Jy to 0.1 Jy, they considered the flat- and 
steep-spectrum populations separately and derived the RLF 
for each population. Their constraints on the RLF of flat- 
spectrum sources were somewhat tighter than those derived 
by Peacock (1985), but the major breakthrough of their pa- 
per concerned the constraints they were able to place on 
the RLF of steep-spectrum radio sources. The strong pos- 
itive evolution in p with z out to z ~ 2 inferred by Lon- 
gair (1966), Wall et al. (1980) and Laing et al. (1983) was 
mapped out in some detail, and more controversially their 
results indicated a decline in the co-moving number density 
of both populations beyond a redshift of m 2 - the so-called 
'redshift cut-off'. DP90 concluded that for powerful sources 
a model of pure luminosity evolution (PLE - evolution only 
of the break luminosity with redshift, analogous to that of 
Boyle, Shanks & Peterson 1988 for the quasar optical lumi- 
nosity function) fitted the data well (but only for Q — 1). A 
similar model which incorporates negative density evolution 
at high-redshift (the LDE model) was found to work well for 
both fi = 1 and Q, — 0. A recent update on the DP90 work 
can be found in Dunlop (1998), and a critical re-evaluation 
of the evidence for a redshift cut-off in the flat-spectrum 
population can be found in Jarvis & Rawlings (2000). 

Derivations of the luminosity functions of luminous 
AGN selected at all frequencies (radio, optical and X-ray) 
have shown a broken power-law form with a steeper slope 
at the high-luminosity end than at the low luminosity end 
(e.g. DP90, Boyle et al. 1988, Page et al. 1996). These studies 
were also able to adequately describe the positive evolution 
of the luminosity function from z — to z m 2 by assuming 
pure luminosity evolution. More recent studies with different 
datasets have shown some deviations from PLE (e.g. Gold- 
schmidt & Miller 1998, Miyaji, Hasinger & Schmidt 2000). 
The physical rationale behind PLE is that the AGN have 
long lifetimes (comparable with the Hubble Time) and they 



decline in luminosity from z ~ 2 to z — 0. However, at least 
in the case of the high-luminosity radio population, we know 
such a model is not viable because the ages of FRII radio 
sources in the surveys are well-constrained to be ^ a few 
10 s yr (e.g. Blunder!, Rawlings & Willott 1999). Further- 
more, for X-ray and optically-selected quasars, accretion at 
the rates required for their high luminosities (over a Hubble 
Time) would lead to much more massive black holes than 
are observed in the local Universe (e.g. Cavaliere & Padovani 
1989). Thus, although PLE models do fit the gross features 
of the data, the physical basis of these models is not clear. 

The radio source population is composed of several 
seemingly different types of object. One division between 
types is concerned with radio structure which falls into two 
distinct classes: FRI, or 'twin-jet', sources; and FRII, or 
'classical double', sources (Fanaroff & Riley 1974). At low 
luminosities log 10 (Lisi) ^25.5, the vast majority of sources 
have an FRI radio structure (apart from the very low lu- 
minosity starbursts), whilst at higher luminosities the ma- 
jority of sources have FRII structure. A second division is 
concerned with the optical properties of the population. At 
log 10 (Li5i) ~ 26.5 the fraction of objects in low-frequency 
selected samples with observed broad-line nuclei changes 
rapidly from w 0.4 at higher radio luminosities to « 0.1 
at lower luminosities (Willott et al. 2000). It is not yet clear 
whether this represents a fundamental change in central en- 
gine properties because, for example, it could be due to the 
opening angle of an obscuring torus increasing with decreas- 
ing luminosity (see Willott et al. 2000). However, it is clear 
that at least some FRII sources lack even indirect evidence 
for an active quasar nucleus because, as pointed out by Hine 
& Longair (1979) and Laing et al. (1994), they have weak or 
absent emission lines. Indeed it is plausible that the major- 
ity of FRII radio galaxies (per unit comoving volume) show 
passive elliptical-galaxy spectra (Rixon, Wall & Benn 1991) 
and are part of the low-luminosity population considered 
here. 

In this paper we have chosen to investigate the RLF 
in terms of a dual-population model where the less radio- 
luminous population is composed of FRIs and FRIIs with 
weak/absent emission lines, and the more radio-luminous 
population of strong-line FRII radio galaxies and quasars. 
Our motivation for doing this is that the presence/absence 
of a quasar nucleus, as indicated by emission line strength, 
is a distinction intimately connected to the properties of 
the central engine, whereas radio structure is strongly influ- 
enced by the larger-scale environment (e.g. Kaiser & Alexan- 
der 1997). The critical radio luminosity [log 10 (-Zjigi) ~ 26.5] 
which crudely divides the populations is also very close to 
the break in the RLF (e.g. from DP90), whereas the luminos- 
ity of the FRI/FRII divide is clearly well below the break. 
Tentative evidence for an increase in the fraction of quasars 
with redshift from z — to 2 w 1 can then be explained if 
this low-luminosity population evolves less rapidly than the 
high-luminosity population (Willott et al. 2000). The p of 
FRI sources is already known to evolve less-strongly with z 
than that of the radio-luminous FRIIs (e.g. Urry & Padovani 
1995), but here we will assume that this is because they 
are a sub-set of a larger population of low-luminosity radio 
sources, including FRIIs, which evolve less strongly. We will 
show that differential evolution between the two populations 
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mimics PLE to a certain extent, but with a more plausible 
physical basis. 

Padovani & Urry (1992), Urry & Padovani (1995) and 
Jackson & Wall (1999) have adopted a different approach to 
determining the RLF from high-frequency complete sam- 
ples. Instead of separating the flat- and steep-spectrum 
sources and deriving the RLF for each population, they use 
models built around unified schemes and relativistic beam- 
ing to determine the relative fractions of various classes of 
objects in the samples. Their models assume that BL Lac 
objects are the beamed versions of FRI and (in the case of 
Jackson & Wall) weak-line FRII radio galaxies, and that 
flat-spectrum quasars are favourably oriented FRII radio 
galaxies. However, both these analyses required the RLF 
of unbeamed objects to be determined independently, i.e 
from low-frequency samples. Because they used only the 
3CRR sample to determine the luminosity distribution of 
radio sources, the tight Li78 — z correlation meant a de- 
generacy between luminosity and redshift effects. Current 
models of this type are thus most useful for constraining 
the values of various beaming parameters and how these af- 
fect the population mix in high-frequency selected samples, 
and have not provided new information on the evolution of 
steep-spectrum sources. 

In this paper we present a new investigation of the 
steep-spectrum RLF which combines the bright 3CRR sam- 
ple with the new, fainter 6CE and 7CRS redshift surveys 
(based on the 6C and 7C radio surveys, respectively). Note 
that an investigation of the RLF of steep-spectrum quasars 
(based on the 3CRR and 7CRS samples) has been presented 
by Willott et al. 1998 (W98). As the basis for an investi- 
gation of the steep-spectrum RLF, the 3CRR/6CE/7CRS 
dataset improves on that used by DP90 in four different 
ways. 

• Spectroscopic completeness. Many of the sources in 
the DP90 samples do not have spectroscopic redshifts [in- 
cluding, at the time of DP90, ~ 50% of the sources in their 
faintest sample - the Parkes Selected Regions (PSR)], and 
redshifts had to be estimated from the infrared Hubble dia- 
gram (the K — z relation). Dunlop (1998) reports that fur- 
ther spectroscopy of these sources has shown their redshifts 
tended to be overestimated (due to the positive correlation 
between radio luminosity and K-banA luminosity, Eales et 
al. 1997), which he contends strengthens the DP90 evidence 
for a redshift cut-off in high-frequency selected samples. 
However, as we discuss later (Section 5.1), the assignment 
of systematically high redshifts to many sources may in fact 
have caused DP90 to overestimate the strength of any red- 
shift cut-off. In contrast to the PSR, the 6CE and 7CRS 
have virtually complete spectroscopic redshifts. 

• Coverage of the radio luminosity versus z plane. 
The flux- limit of the PSR is S2.7 = 0.1 Jy. For steep- 
spectrum sources with a = 0.8, this corresponds to £151 ~ 1 
Jy, and for a = 1.2toS'i5i~3 Jy, so the 7C Redshift Sur- 
vey (7CRS) with a flux-limit of S151 = 0.5 Jy is a factor 
2 — 6 times deeper at a given z than the PSR depending on 
spectral index. Because of the correlations between a and 
radio luminosity and z (e.g. Blundell et al. 1999) this makes 
the 7CRS sensitive to lower radio luminosity objects at the 
higher redshifts. Note that because the PSR covers an area 



3.5 times larger than the 7CRS, the two samples actually 
contain similar numbers of high-redshift objects. 

• Radio completeness. DP90 found the source counts 
in the PSR at 0.1 Jy to be a factor of two lower than expected 
from the source counts at other frequencies. They suggested 
that this is due to incompleteness in the original survey near 
its flux-limit, although this is by no means certain (J.V. 
Wall, priv. comm.). Although they attempted to correct for 
this, there remains some possibility of bias: methods based 
on the y/Vmax statistic are particularly sensitive to such 
problems, as objects close to the flux- limit have V/Vma* ~ 1, 
and losing such objects would bias the statistic towards a 
spuriously low value. The 6CE sample and the 7CRS are 
based on surveys which are complete at the flux-densities 
of interest. As Riley (1989) discusses, the number of large 
angular size sources omitted from surveys such as these due 
to surface brightness effects is likely to be negligible. 

• Benefits of low— frequency selection. The ultimate 
goal of our programme is to model the evolution in the RLF 
of steep-spectrum sources. As detailed in Blundell et al. 
(1999), the physics of FRII radio sources means that the 
mapping between radio luminosity and one of the key phys- 
ical variables, the bulk power in the jets, is much closer if 
the radio luminosity is evaluated at low rest-frame frequen- 
cies provided, of course, these lie above the synchrotron self- 
absorption frequency of a given source. Modelling of high 
frequency luminosity is much more problematic, both be- 
cause lobe emission at these frequencies is highly sensitive 
to source age and redshift (because of inverse-Compton scat- 
tering of lobe electrons by the microwave background) , and 
because Doppler-boosted cores and hotspots become much 
more important at high rest-frame frequencies. 

In Section 2 we discuss the data available to constrain 
the steep-spectrum RLF. Section 3 describes the model fit- 
ting procedures and results. In Section 4 the evolution of the 
radio source population is tested via the V/Knax statistic. 
Section 5 compares our findings with those of previous RLF 
determinations. In Section 6 we summarise our results and 
indicate the future work necessary to resolve the outstanding 
uncertainties. 

The convention for the radio spectral index, a rat j, is 
that S v oc i>~ a ™ d , where Si, is the flux-density at frequency 
v. We assume throughout that H = 50 km s -1 Mpc~ and 
Qa = 0. All results are presented for both flat (Qm = 1) and 
open (^m = 0) cosmologies. 



2 CONSTRAINING DATA 

The radio luminosity function (RLF) is defined as the num- 
ber of radio sources per unit co-moving volume per unit 
(base 10) logarithm of luminosity, p(L,z). To determine the 
RLF, we will make use of three different types of radio 
data: complete samples with spectroscopic redshifts; source 
counts; and a low-redshift sample to fix the local RLF. 



2.1 Complete samples 

A complete sample consists of every radio source in a cer- 
tain area on the sky brighter than a specified flux-limit at 
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Figure 1. The radio luminosity-redshift plane for the 3CRR, 6CE and 7CRS samples described in Section 2.1. The different symbols 



identify sources from different samples: 3CRR (filled circles), 6CE (open diamonds) and 7CRS (asterisks). The dotted lines show the 
flux-density lower limits for the 3CRR and 7CRS samples and the dashed lines the limits for the 6CE sample (all assuming a radio 
spectral index of 0.5). This plot is for (Im = 1, = 0. 



the selection frequency. If all these sources are identified op- 
tically with a galaxy or quasar and redshifts are determined, 
then the distribution of points on the L — z plane, together 
with knowledge of the sample fiux-limits and sky areas, in 
principle allow the determination of the RLF. The major 
uncertainty in such a determination is that at least some 
regions of the L — z plane will be very sparsely populated, 
and at these L and z the RLF measurement will be prone to 
small number statistical uncertainties. These typically arise 
because the identification of and/or redshift acquisition of 
faint radio sources can be very time-consuming, and due to 
the limited nature of time available at large telescopes, com- 
plete samples up until now have generally been limited to 
;$ 100 sources. Some regions of the L — z plane are prone 
to a more fundamental problem: objects of a given L may 
be sufficiently rare that there is simply too little observable 
comoving volume at the redshift of interest to obtain large 
samples. 

Our dataset of redshift surveys combines three low- 
frequency complete samples with different flux-limits, which 
together provide unparalleled coverage of the L151 — z plane 
(see Fig.M. 



2.1.1 ICRS 

The faintest of the three complete samples used is the new 
7C Redshift Survey (7CRS). This sample contains every 
source with a low-frequency flux-density Si 51 > 0.5 Jy in 
three regions of sky covering a total of 0.022 sr. The regions 
7C-I and 7C-II will be described in Blundell et al. (in prep.) 
and 7C-III in Lacy et al. (1999a) and references therein. For 
7C-I and 7C-II, multi-frequency radio data including high- 
resolution radio maps from the VLA will be presented in 
Blundell et al. From measurements at several frequencies, 
the radio spectra have been fitted and hence radio luminosi- 
ties at rest-frame 151 MHz determined. Spectral indices at 
various frequencies were also determined from these fits. For 
7C-III, the availability of 38 MHz data from the 8C survey 
in the same region has allowed the low-frequency spectral 
indices to be calculated and hence rest-frame 151 MHz lu- 
minosities (Lacy et al. 1999a). 

Due to the low-frequency selection of the sample, it 
is expected to contain few flat-spectrum (chghz < 0.5) 
sources. From inspection of the VLA maps and radio spec- 
tra, sources which only have flux-densities above the sam- 
ple limit due to their Doppler-boosted cores can be iden- 
tified. Out of the total of 130 sources, only one source 
fulfills these criteria - the quasar 5C7.230. Therefore this 
source is excluded from the analysis presented in this paper. 
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Note however that there are five other quasars which have 
aiGHz < 0.5, the traditional separation value between flat- 
and steep-spectrum quasars. All of these have qighz ~ 0.4 
and are compact (6 < 5 arcsec), but have extended emission 
on both sides of the core, indicating that Doppler-boosting 
is not the cause of this emission. They all have sufficient 
extended flux to remain in the sample and as such they are 
not excluded on orientation bias grounds. These sources are 
referred to as Core-JetS sources (CJSs; see Blundell et al. 
in prep, for more details). The one other exclusion is 3C200, 
because it is already included in the 3CRR sample, leaving 
a total of 128 sources. 

All 7C Redshift Survey sources have been reliably iden- 
tified with an optical/near-infrared counterpart. For 7C-I 
and 7C-II, spectroscopy has been attempted on all sources 
and secure redshifts obtained for > 90% of the sample 
(Willott et al., in prep.). Seven sources do not show any 
emission or absorption lines in their spectra. For these 
objects, we have obtained multi-colour optical and near- 
infrared photometry (in the R, I, J, H and K bands) to 
attempt to constrain their redshifts. All of these objects have 
spectral energy distributions consistent with evolved stellar 
populations at redshifts in the range 1 < z < 2 (Willott, 
Rawlings & Blundell 2000) . These derived redshifts are gen- 
erally consistent with those from the K — z diagram, partic- 
ularly given the increased scatter in the relationship at faint 
magnitudes (Eales et al. 1997). Note that radio galaxies in 
the redshift range 1.3 < z < 1.8 are traditionally difficult 
to obtain redshifts for because there are no bright emission 
lines in the optical wavelength region. The correlation be- 
tween emission line and radio luminosities (Baum & Heck- 
man 1989; Rawlings & Saunders 1991; Willott et al. 1999) 
exacerbates this problem for faint radio samples. We believe 
that this is the primary reason why we have been unable 
to detect emission lines in these 7 sources. It is extremely 
unlikely that any of these objects actually lie at z > 2.5. 

Details of the imaging and spectroscopy in the 7C-III 
region are given in Lacy et al. (1999a, 1999b, 2000). Rea- 
sonably secure spectroscopic redshifts have been obtained 
for 81% of the sources and uncertain redshifts for a further 
13%. Near-infrared magnitudes recently obtained support 
these uncertain redshifts. Only three sources have no red- 
shift information. Near-infrared and optical imaging shows 
that two of them have SEDs typical of z ~ 1.5 galaxies as 
is the case for those without spectroscopic redshifts in 7C-I 
and 7C-II. The other object (7C 1748+6703) has a break be- 
tween J and K that is suggestive of a redshift z > 2.4. The 
existence of significant flux in its optical spectrum down to 
6000 A and a isT-magnitude of 18.9, suggest z < 4. Therefore 
we adopt z — 3.2 in this paper. 

2.1.2 3CRR 

The bright complete sample used is the 3CRR sample of 
Laing, Riley & Longair (1983, hereafter LRL) with complete 
redshift information for all 173 radio sources. The flux-limit 
is 10.9 Jy at 178 MHz, which translates to 12.4 Jy at 151 
MHz assuming a typical spectral index of 0.8. Blundell et 
al. (in prep.) review the current status of this sample and 
presents the results of spectral fitting of multi-frequency ra- 
dio data. 3C231 (M82) is excluded because the radio emis- 
sion from this source is due to a starburst and not an AGN. 



The two flat-spectrum quasars, 3C345 and 3C454.3, are ex- 
cluded because only Doppler-boosting of their core fluxes 
raises their total fluxes above the flux-limit. To include ob- 
jects such as these would cause an overestimate in the num- 
ber density of AGN, since they only get into the samples 
because their jet axes are very close to our line-of-sight. 
Note that high-frequency selected samples contain a much 
higher proportion of objects such as these, showing that 
our choice of low-frequency selection virtually eliminates 
these favourably oriented sources [there are only three such 
sources (~ 1%) in our complete samples]. 



2.1.3 6CE 

The 6CE complete sample we use is a revision of the 6C 
sample of Eales (1985) (see Rawlings, Eales & Lacy 2000 
for details). The flux-limits of this sample are 2.0 < S151 < 
3.93 Jy and the sky area covered is 0.103 sr. Only 3 of the 
59 sources in this sample do not have redshifts determined 
from spectroscopy. One object is occluded by a bright star 
and is excluded from further consideration (without bias). 
One object is faint in the near-infrared (K > 19), so we 
take it as a galaxy at z = 2.0 in this paper. The other is 
relatively bright in the near-IR but very faint in the optical. 
Its red colour suggests a galaxy with a redshift in the range 
0.8 < z < 2 and we assume z — 1.4 here. None of the 6CE 
sources are clearly promoted into the sample by Doppler- 
boosted core emission (c.f 5C7.230, 3C345 and 3C454.3), so 
all 58 sources are included in our analysis. Full details of the 
sample, including optical spectra, are given in Rawlings et 
al. (2000). 



2.2 The source counts at 151 MHz 

The three complete samples with known redshift distribu- 
tions contain a total of 356 sources (after the exclusions 
mentioned above). However, the number of radio sources at 
151 MHz as a function of flux density, also known as the 
source counts, has been determined for thousands of sources 
from the much larger sky area 6C and 7C surveys (Hales, 
Baldwin & Warner 1988 and McGilchrist et al. 1990, respec- 
tively). The 7C source counts go as faint as 0.1 Jy - a factor 
of five fainter than the 7CRS. At these flux-densities the 
source counts are our only constraint on the RLF. At the 
bright end ( > 10 Jy) , the source counts are obtained by bin- 
ning the sources in the 3CRR sample of L RL in 5 flux bins. 
The three sources excluded in Section 2.1 were also excluded 

The 3CRR source 



for the calculation of the source counts 
counts are converted from 178 to 151 MHz by assuming a 
spectral index of 0.8. 

Figure ^ shows the binned differential source counts 
from the 3CRR, 6C and 7C surveys. The counts are nor- 
malised to the differential source counts for a uniform 
distribution in a Euclidean universe, such that dNo = 
2400 (S~l n 5 - S~lx), where S m in and S max are the lower 
and upper flux limits of the bin. These points were fitted 
by a third-order polynomial using a least-squares fit in log- 
arithm space, i.e. 



fdN\ 
l0Sl0 [dNo ) 



a + en log 10 S151 + 

02 (log 10 S151) 2 + a 3 (log 10 S151) 3 . 



(1) 
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Figure 2. The binned differential source counts from the 3CRR 
(LRL), 6C (Hales et al. 1988) and 7C (McGilchrist et al. 1990) 
surveys. Error bars are Poisson errors for the number of 

sources in each bin. The symbols identify the surveys: 3CRR (cir- 
cles), 6C (crosses) and 7C (triangles). The line shows the third- 
order polynomial fit to the data points. 



ties below L 151 « 10 24 W Hz -1 sr~ x . Hence the LRLF will 
only constrain the faint end of the RLF. Another problem 
is that at very low luminosities (L151 « 10 21 W Hz _1 sr _1 ) 
the LRLF is dominated by the star-forming spiral/irregular 
galaxies which dominate the sub-mjy (at 1.4 GHz) source 
counts. Since we do not want to include this population in 
our RLF determination, it is essential that the LRLF of the 
AGN population only is used. 

A determination of the local radio luminosity function 
which treats AGN and star-forming galaxies as two distinct 
populations is given by Cotton & Condon (1998). By cross- 
matching UGC galaxies (Nilson 1973) with the NVSS cat- 
alogue (Condon et al. 1998) they derived luminosity func- 
tions at 1.4 GHz for each population. The AGN LRLF at 
151 MHz was derived from this by assuming a radio spec- 
tral index of 0.8. This determination of the LRLF by Cot- 
ton & Condon is well-approximated by a power-law (up to 
L151 « 10 24 W Hz" 1 sr" 1 ) with an index of -0.53. From 
this model the LRLF was determined in 10 bins in the range 
20 < log 10 (Li5i/ W Hz" 1 sr" 1 ) < 23.6. The error on each 
bin was assumed to be 0.1 dex. 



The coefficients of the fit are ao = —0.00541, a\ — 
0.293, a 2 = -0.362 and a 3 = 0.0527. A higher order poly- 
nomial was not required to provide a decent fit to the data. 
The error on this curve was estimated using the errors from 
the source count bins. At 6*151 < 1.0 Jy the fractional errors 
on the bins are approximately constant with flux-density at 
the level of 0.06. At S151 > 1.0 Jy the fractional errors are 
given by a power-law with slope 0.49 such that the fractional 
error at 10 Jy is 0.19, for example. 

It is well known that there is a flattening of the source 
counts slope at 1.4-GHz flux-densities ~ 1 mjy. This is at- 
tributed to the counts becoming dominated by low-redshift 
star- forming galaxies (Condon 1984; Windhorst et al. 1985). 
Although there is only one such starburst galaxy in our com- 
plete samples down to 0.5 Jy (3C231 ), the source count data 
goes as low as 0.1 Jy, so it is essential to investigate the con- 
tribution of starburst galaxies to the source counts at these 
low flux-densities. To estimate the starburst contribution we 
used the 60/xm luminosity function of Saunders et al. (1990) 
scaled by the mean of the ratio of IRAS 60/im to NVSS 1.4 
GHz flux-densities (Cotton & Condon 1998). Cosmic evolu- 
tion is accounted for by the PLE model adopted by Saunders 
ct al. To translate this 1.4 GHz luminosity function to 151 
MHz, a spectral index of 0.7 was assumed for the starburst 
population. The differential source counts at 151 MHz were 
then calculated from this luminosity function. It is found 
that the starburst contribution at 0.1 Jy is only 2% of the 
total source counts. This is within the uncertainty in the 
source counts at 0.1 Jy, so it can be safely neglected. 



2.3 The local RLF 

The radio luminosity function at low-redshift (z < 0.2) is 
fairly well-determined. This therefore provides a boundary 
condition for determination of the evolving RLF. Unfortu- 
nately, the low space density of powerful radio sources and 
small cosmological volume out to z w 0.2 means that the lo- 
cal RLF (LRLF) is only well-constrained at radio luminosi- 



3 MODEL FITTING 
3.1 Method 

To find best-fit parameters for various models of the evolv- 
ing RLF, a maximum likelihood method is adopted. This 
method is similar to that used in W98 to determine the 
quasar RLF and we refer readers to that paper for further 
details. The aim of the maximum likelihood method is to 
minimise the value of S, which is defined as 



S= -2^2\n[p(Li,Zi)] + 



'If 



p(L, z)fl(L, z)^dzdlog 10 L, 



(2) 



where (dV/dz)dz is the differential co-moving volume ele- 
ment, fl(L, z) is the sky area available from the samples for 
these values of L and z and p(L, z) is the model distribution 
being tested. In the first term of this equation, the sum is 
over all the N sources in the combined sample. The second 
term is simply the integrand of the model being tested and 
should give w 2N for good fits. 

The only difference here is that the source counts and lo- 
cal radio luminosity function provide additional constraints 
for the fitting process. The source counts and LRLF are one- 
dimensional functions and to estimate how well they fit the 
model, the value of x 2 is evaluated thus: 



N 
i=l 



/"data i fmod 
^"data i 



■\ 2 

1 ■ 



(3) 



where /data ; is the value of the data in the ith bin, similarly 
/mod i and Odata i are the model value and data error in the 
ith bin, respectively. These are all determined in logarithm 
space for both the source counts and the LRLF. The sum is 
over the N bins in each case. 

In order to combine the constraints from all three types 
of data, we follow Kochanek (1996). \ 2 ls related to the 
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likelihood by \ 2 = —2 In (likelihood), i.e. the same form as 
S. Therefore the function that is now minimised is 



Sail — Xsc + Xlrlf + S 



•So, 



(4) 



where So is a constant which normalises S so that equal sta- 
tistical weight is given to all three types of data. The value 
of So theoretically comes from the dropped terms from equa- 
tion [| However, these are not possible to calculate directly, 
so instead one estimates the value of So so that the first 
term in equation ^ equals « and S w 2iV. Errors on the 
best-fit parameters are 68% confidence levels, as discussed 
in Boyle et al. (1988). 

The goodness-of-fit of models are estimated using the 
2-D Kolmogorov-Smirnov (KS) test (Peacock 1983), as in 
W98. Note that it is expected that the KS probabilities will 
be lower than for the quasar RLF fits, because in general 
the probability decreases with increased sample size in this 
test. Simulations of the L151 — z plane are used to double 
check the goodness-of-fit. In addition, the predicted redshift 
distributions for the complete samples are derived from the 
models and compared with the actual distribution. This is 
particularly useful in testing for the presence of a redshift 
cut-off at high-z. 



3.2 The form of the RLF 

As explained in Sec. [j], the RLF is modelled here as a com- 
bination of two populations which are allowed to have dif- 
ferent RLF shapes and evolutionary properties. The low- 
luminosity radio source population (composed of FRIs and 
low-excitation/weak emission line FRIIs) is modelled with 
a Schechter function, 



pi(L) = pio 



(5) 



where p\(L) is the source number density as a function of 
radio luminosity L, p\ is a normalisation term, Lu is the 
break luminosity and ol\ is the power-law slope. The evolu- 
tion function of p\ is modelled as fi(z) — (1 + z) kl up to 
a maximum redshift z\ beyond which there is no further 
evolution. Hence there are five free parameters to be fixed 
for the low-luminosity population. 

For the high-luminosity population a similar form is 
adopted, except in this case the exponential part of the 
Schechter function is at low luminosity and the power-law 
region is at high luminosity, i.e. 



ph(L) = p hi 



L 



(6) 



where the subscript h refers to the high luminosity popu- 
lation. Note that it is expected that Li* ~ Lh*, so that 
at this luminosity the decline in one population is compen- 
sated by the rise of the other, since no discontinuity or bi- 
modality in the RLF has yet been observed (e.g. DP90). 
These forms for the luminosity functions of the two popu- 
lations were chosen because the steep exponentials ensured 
a fairly small overlap in L151 between the two populations. 
However, we currently have little information about the rel- 
ative population fractions in the overlap region, so this is 
just a convenient method of separating the populations. We 
note that the low-luminosity RLF of our model has a sim- 



ilar shape to the elliptical galaxy _K"-band luminosity (and 
therefore mass) function (e.g. Gardner et al. 1997). 

The high-luminosity population is assumed to undergo 
density evolution of a similar form to that of the quasar RLF 
in W98. The z-distribution fh(z) is given by a Gaussian in 
redshift (model A here), with a variant on this distribu- 
tion being a one-tailed Gaussian from zero redshift to the 
peak redshift and then a constant density to high-redshift 
(model B - also known as the 'no cut-off' model since there 
is no redshift cut-off). Note that the Gaussian model A is 
a convenient approximation to the shape of the evolution of 
luminous optically-selected quasars (e.g. Warren, Hewett & 
Osmer 1994) , and to the evolution inferred for flat-spectrum 
radio sources by Shaver et al. (1999) - but see Jarvis & 
Rawlings (2000) for a contrary view. There is no physical 
reason why one would expect the evolution to be symmetric 
in redshift. In W98 we showed that the Gaussian redshift 
distribution from z — to the peak is very similar to the 
form (1 + z) k , so we do not include models with this form 
of evolution here. 

In W98 it was found that both these redshift distribu- 
tions (models A and B in this paper) provided good fits be- 
cause the number of quasars in the samples at high-redshift 
was fairly small. The inclusion of the 6CE and 7C-III sam- 
ples and the 7CRS radio galaxies gives a much larger sample 
of 2 > 2 objects here than for the quasar RLF. Hence, the 
high redshift evolution can be better determined. Therefore, 
an additional redshift distribution is tested here, which has 
a one-tailed Gaussian rise to the peak redshift and then a 
one-tailed Gaussian decline at higher redshifts which is al- 
lowed to have a different width from the rise. This allows 
more shallow high-redshift declines in ph (such as that in- 
ferred by DP90) than the symmetric model A. This model 
(C) directly fits the strength of any decline in the co-moving 
density beyond the peak redshift Zho, since the value of the 
high-redshift Gaussian width z\a is free to be fit from for 
an abrupt cut-off to a large value for a very gradual decline 
at high-redshift. Hence there are a total of six free parame- 
ters in the high-luminosity population model with this form 
of redshift evolution. 

To summarise our adopted RLF models, below we give 
all the equations needed to reproduce the RLF p(L, z). The 
values of the best-fit parameters are given in Table ^. 



p(L, z) = pi + p h 
where 

L \ " ai ( —L 



(7) 



Pi = Pi 



Pi 



o( — ) ' cx p( — )( 1 + z )" 1 forz<zio, (8) 
Pio ^ — J exp ^ — J (1 + zio) 1 for z > z\ , (9) 



The high-luminosity evolution function /h(z) has three dif- 
ferent forms depending upon the model and the redshift: 



/ h (z) = exp 



1 / z - z hc 



2 V Zhi 



(11) 



for model A at all redshifts and models B and C at z < z^ ; 
Mz) = 1.0 (12) 
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for model B at z > Zho 



/h(z) = exp 



1 / z - Zho y 

2 I Z h2 / 



(13) 



for model C at z > Zho. 

In this paper, we have performed the RLF modelling for 
two cosmological models with Qm = 1 and Qm = and zero 
cosmological constant. To obtain a fairly reliable estimate 
of the RLF in other cosmologies one can use the following 
relation from Peacock (1985), 



it \ dVl n \ dV2 
Pi{Li,z)— = p 2 (L 2 ,z) — , 

az az 



(14) 



where L\ and L2 are the luminosities derived from the flux- 
density and redshift in the two cosmologies. We have tested 
this relation by converting our model C RLF derived for 
Qm = back to Qm = 1 and then compared this with 
the RLF derived for Qm = 1. We note that the relation 
works well for regions of the luminosity function that are 
well-constrained by data. However, in regions which are con- 
strained by little or no data, the relation may not reproduce 
a model RLF close to that which would emerge from di- 
rect fitting to the data in the new cosmology. Thus, Table [I] 
shows that the model C peak redshifts Zho and strength of 
high-z decline Zh2 are different for Qm = 1 and Qm = and 
equation |l4| cannot reproduce such differences. Note that 
the RLF for Qm = 0.3, Qa = 0.7 is very similar to that of 
Qm ~ 0, Qa ~ for regions of the RLF well-constrained by 
our data and equation ^ can be applied in this case. 



3.3 Results 

The maximum likelihood fitting routine was run for these 
three RLF models. The routine integrated over the ranges 
20 < log 10 (I/151 ) < 30 and < z < 5. A check was made to 
ensure that the same results occurred if the redshift range 
was increased to z = 10. All three models provided good 
fits to the data. The best-fit parameters for each model are 
shown in Table [jj 

The parameters for the low-luminosity population RLF 
do not change significantly for the different models, which 
is as expected, since it has little to do with the form of 
the high-z evolution of the high-luminosity RLF. Note that 
the value of the break luminosity, log(Li*) ~ 26.1 (26.5 for 
Qm — 0), is above the FRI/FRII divide and below the break 
in the quasar RLF from W98 [log(L brcak ) w 26.8 (27.2 for 
Qm = 0)]. The models show significant evolution (fei ta 4) 
for the low-luminosity objects out to z ~ 0.7 with no further 
evolution to higher redshifts. Although it is often stated that 
the FRI radio galaxy population shows no evolution, the 
samples used to determine this generally have small numbers 
of objects and a narrow redshift range (e.g. Urry, Padovani 
& Stickel 1991; Jackson & Wall 1999). As we shall see in 
Section | the combination of the 3CRR, 6CE and 7CRS 
samples shows considerable low-redshift evolution for low- 
luminosity (log 10 ii5i < 26.5) FRI and FRII radio galaxies. 

The best-fit models for the high-luminosity RLF have 
much steeper slopes than the low-luminosity RLF (a?h ~ 2.3, 
c.f. Qi « 0.6). This, together with log(Li*) ~ log(Lh*), gives 
the whole RLF a broken power-law form, such as observed 
by DP90 for the RLF and by Boyle et al. (1988) for the 



Model 


Qm 


iV par 


Pks 


A 


1 


10 


0.11 


13 


1 


10 


0.28 


C 


1 


11 


0.36 


A 





10 


0.03 


13 





10 


0.09 


C 





11 


0.07 



Table 2. Goodness-of-fit table. -/V par is the number of free pa- 
rameters in the model (including those from both the low- and 
high-luminosity populations) and Pks is the probability from the 
2D K-S test to the L151 — z distribution. As in Table [l] the top 
three rows refer to Qm = 1 and the bottom three to Om = 0. 



quasar OLF. Lh* in these models is about an order of mag- 
nitude greater than L\+ and is very similar to the break in 
the quasar RLF estimated in W98 from source counts con- 
straints. The RLF of model C at various redshifts is shown in 
Figure ^| The broken power-law form for the whole popula- 
tion can clearly be seen at z < 1. At higher redshifts the con- 
tinued evolution of the high-luminosity population causes 
the RLF to have a small dip at log(Lisi) ~ 26.2 (26.6 for 
Qm = 0). This is a consequence of the models chosen and 
may or may not be real. No previous low-frequency com- 
plete samples have gone to such low flux-density levels to 
constrain this region of the L151 — z plane directly. In fact, 
for log(Lisi) ~ 26.2 and z > 1, the 151 MHz flux-density 
of any sources would be fainter than the limit of the 7CRS 
sample and the RLF is purely constrained by the source 
counts in this region. The maximum size of this feature is 
only 0.2 dex, so it would be extremely difficult to confirm 
directly by measurement of the RLF. 

Fig. |j shows the calculated source counts at 151 MHz 
from the model C RLFs for Qm = 1 and Qm = 0. The counts 
from each of the two parts of the RLF have also been plotted 
separately. The models fit the total source counts very well, 
as is expected because the maximum likelihood method min- 
imises the x 2 °f the fit to the source counts data. For both 
cosmologies, the high-luminosity population dominates the 
source counts over most of the flux-density range. The low- 
luminosity population starts to dominate at ~ 0.2 Jy. This 
suggests that fainter samples, such as those selected at flux- 
limits 5*151 > 0.1 Jy, would contain more low-luminosity 
population sources. This idea will be discussed further in 
Section |(| 

We now briefly consider whether it is obvious that any 
of these three models provides the best-fit to the data and 
hence whether we can usefully constrain the strength of any 
high-z decline in the co-moving space density p of radio 
sources. We defer a full statistical analysis of this question 
to a future paper (Jarvis et ah, 2000) in which data from 
the 3CRR and 6CE complete samples is combined with the 
results of targeted searches for z > 4 radio galaxies (the 
6C* sample of Blundell et al. 1998). Figure M shows the 
redshift-dependence of the high-luminosity RLF for these 
three models with their free parameters fixed at the best-fit 
values. All the models provide good fits to both the source 
counts and local RLF data, with reduced \ 2 values of order 
one. Table ^| shows the goodness-of-fit (KS test probabilities) 
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Figure 3. The radio luminosity function derived for model C for !)m = 1 (top panel) and Qm = (bottom panel). Dashed lines show the 
low-luminosity population at z = 0, z = 0.5 and z ?J 0.7 (bottom to top). The dotted lines are the high-luminosity RLF at z = 0, 0.5, 1, 3 
and 2 (bottom to top). Solid lines show the sum of both components. 
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Model 


Qm 


log(pio) 


ay 


log(ii*) 


Zlo 


fci 1 


Og (Pho 


) cth 


log(L h *) 


Zho 




2h2 


A 


1 


-7.153 


0.542 


26.12 


0.720 


4.56 


-6.169 


2.30 


27.01 


2.25 


0.673 


- 


B 


1 


-7.150 


0.542 


26.14 


0.646 


4.10 


-6.260 


2.31 


26.98 


1.81 


0.523 


- 


C 


1 


-7.120 


0.539 


26.10 


0.706 


4.30 


-6.196 


2.27 


26.95 


1.91 


0.559 


1.378 


A 





-7.503 


0.584 


26.46 


0.710 


3.60 


-6.740 


2.42 


27.42 


2.23 


0.642 




B 





-7.484 


0.581 


26.47 


0.580 


3.11 


-6.816 


2.40 


27.36 


1.77 


0.483 




C 





-7.523 


0.586 


26.48 


0.710 


3.48 


-6.757 


2.42 


27.39 


2.03 


0.568 


0.956 


Errors for model C (Hm = 


1): log(p 


lo); -0.11, +0.10, 


an -0.02, 


+0.02, 


logfa*); - 


-0.09, +0.08, z lo ; - 


-0.10,+ 


0.10, fci; 


.55, +0.57, 


log (Pho); -0.11, +0.09, 


a h ; -0.11,- 


f0.12, 


log(L h *); 


-0.10,4 


-0.11, Z ho 


-0.16,+ 


0.16, z hl 


; -0.05 


+0.05, z h2 



-0.28, +0.52. 



Table 1. Best-fit parameters for RLF models A, B and C, as described in Section 3.2, The top three rows refer to Om = 1 and the 
bottom three rows Cm = 0. 
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Figure 4. The 151 MHz source counts data (as in Fig. |^) with 
model-B (short-dash) and model-C (long-dash) fits for SJm = 1 
(top) and Om = (bottom). The contributions to the source 
counts from the low-luminosity and the high-luminosity popula- 
tions are shown separately (and marked 'hi' and 'lo'). The solid 
lines show the total source counts from both models. 



of the models for the two cosmologies considered. In both 
cosmologies the KS test shows that models B and C are 
marginally favoured over model A - the symmetric decline 
model. The extra parameter in model C (zh2 _ the slope of 
the decline in number density at high-redshift) enables the 
strength of any high-z decline to be fitted directly. The fact 
that the best-fitting values of Zh2 are greater than those of 
Zhi but less than infinity indicates a moderate high-z decline 



Figure 5. The best-fit redshift dependence of the high- 
luminosity part of the RLF for the 3 models considered here 
(fijyi = 1 - solid lines, S!m = - dashed lines). Note that the 
curves for model B in the two cosmologies are the same at high- 
redshift, so the dashed curve is obscured by the solid curve. All 
are normalised such that log 10 (p) = at the peak redshift. 



fits the data best (see Fig. |B|). We conclude, however, that 
we cannot rule out with any level of confidence a constant 
space density up to z ~ 5. 



3.4 Simulated data 

Simulations of the L151 — z plane for the three models are 
shown in Figure |^ along with the actual L151 — z data for 
the complete samples (for Qm = 1). Statements made in 
this section concerning comparison of real and simulated 
data apply also to an £Im = cosmology. The first thing 
to notice is that all the models reproduce the L151 — z 
data reasonably well. Sources with emission line luminosi- 
ties log 10 (L [O ii]/W) > 35.1 (see Willott et al. 1999, 2000 for 
discussion of this adopted value) are shown as filled symbols 
whereas those with lower luminosities are open. Note that 
because the radio galaxies in the combined sample span a 
large range of redshift, Balmer lines are often not observed 
and a classification scheme such as that adopted by Laing et 
al. (1994) to distinguish between low- and high-ionisation 
sources cannot be used. Therefore we make a distinction be- 
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Figure 6. The top-left panel shows the L151 — z plane for the 3CRR, 6CE and 7CRS complete samples. Objects in the samples with 
emission line luminosities log 1 o(^[on]/W) > 35.1 are shown as filled symbols whilst those at lower luminosities are open. The other 
three panels show simulations of the L151 — z plane for samples with the same flux-limits and sky areas as the complete samples. These 
simulations were generated using the best-fit RLF models A, B and C (S!m = 1). Objects drawn from the low— luminosity population are 
shown as open circles and those from the high-luminosity population are open. 



tween the two populations simply on the basis of emission 
line luminosity. 

For the model simulations, the high and low luminos- 
ity populations are similarly separated into filled and open 
symbols. The transition between populations at log 10 (Lisi/ 
W Hz -1 sr _1 ) as 26.5 is well-reproduced by the models indi- 
cating that the two populations can approximately be sep- 
arated according to the strength of their narrow emission 
lines. As seen in Willott et al. (1999), this quantity is likely 
to be closely related to the accretion rate of the source. The 
scatter in the emission line-radio correlation is not taken 
account of in the RLF models, which explains why the filled 
and open symbols overlap in the real data, but do not for the 
simulations. This scatter is due to effects such as the range 
of radio source ages and environments which give different 
radio luminosities for a given jet power (see Willott et al. 
1999; Blundell et al. 1999). 

At high-redshifts (2«3) the very few numbers of ob- 
jects in the samples means that small number statistics be- 
come important. For a comparison with simulated data, the 
situation is even worse because both the data and the simu- 
lations have independent Poisson errors. Therefore, the sim- 
ulations shown in this section should not be used to dis- 



tinguish between best-fitting models (a direc t in tegration of 
the RLF over redshift is preferable, see Sec. 3.5), but these 



simulations could give an indication of a poor fit. In any 
case, the simulations look very similar, even at z fa 3, for all 
the three models tested. 



3.5 Redshift distributions 

We now consider the redshift distributions for the 3CRR, 
6CE and 7CRS complete samples. As these are obtained by 
integrating the RLF, the only source of small number statis- 
tics come from the data they are compared with. Figure 
shows histograms of the number of sources per bin of width 
Sz = 0.25 for the 3CRR, 6CE and 7CRS samples separated 
into two populations depending upon emission line luminosi- 
ties as before. For Qm = 0, the high/low luminosity divide 
is at log 10 (L[oii]/W) = 35.4. Also plotted are the model C 
RLF predictions for each population for the three samples. 

For the 3CRR sample, the model fits very well, the main 
exception being the low luminosity objects at z > 0.3, where 
there are several more sources observed than predicted. This 
was also seen in the simulations and attributed to scatter in 
the emission line-radio correlation. With only 58 sources in 
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Figure 7. The histograms show the number of sources in the 
3CRR, 6CE and 7CRS complete samples binned in redshift with 
bin width Sz = 0.25 for Cm = 1 and f!m = 0. There are separate 
histograms for the low (dashed) and high (dotted) luminosity 
populations. The solid lines show the model C predictions for the 
number-redshift distribution for each population. Note that the 
3CRR plot only has a log scale on the y-axis because of the high 
peak at low-redshift. 



the 6CE sample, the effects of small number statistics are 
clear with an irregular observed distribution. However, in 
general the model redshift distribution is a fairly good ap- 
proximation to the data. There is something of an excess of 
low-luminosity sources predicted at low-redshifts [z < 0.5), 
but this can be accounted for by the scatter, since there are 
several high-luminosity objects observed at these redshifts. 
For the 7CRS sample, again we find a fairly good fit. The 
main differences here are a small deficit predicted at z ~ 0.5 
and an excess at z « 1.5. The redshift distributions at high- 
redshift (z > 2) will be discussed later in this section. 

Figure ^ plots the cumulative redshift distributions of 
the two samples for all three models. The similarity of these 
three curves shows that the models all predict similar red- 
shift distributions of sources. The one-dimensional KS test 
was used to test the difference between the observed and 
model redshift distributions. For the 3CRR sample the sig- 
nificances are 0.011 (0.03), 0.007 (0.03) and 0.007 (0.02) for 
models A, B and C, respectively (Qm = values in brack- 
ets). The values for the 6CE sample are 0.40 (0.32), 0.15 
(0.10), 0.16 (0.18) and for the 7CRS sample 0.09 (0.002), 
0.04 (0.0004) and 0.09 (0.001). The lowest values of KS sig- 



Figure 8. Normalised cumulative redshift distributions of 
sources in the 3CRR, 6CE and 7CRS complete samples. The solid 
curves are the cumulative redshift distributions in the samples 
predicted by the RLF of models A (dotted), B (dashed) and C 
(dot-dashed). 



nificance are for the 7CRS redshift distribution for the case 
of J1m = 0. It is fairly easy to speculate on the cause of this 
low probability: the model RLFs appear to be too high at 
Li5i ~ 26.5 and too low at z ~ 0.5 and 



z ~ 1.5 and log 10 liVo\ 



login il5 



26 (fin 



1, for Q,\ 



these luminosities 



are log 10 -Lisi ~ 26.9 and 26.4). This is almost certainly a 
limitation of our modelling procedure since these luminosi- 
ties lie close to the break luminosities Li* and Lh* which 
bracket the range over which the two populations overlap. 
Away from these break luminosities the rates of exponential 
decline in luminosity of each population are fixed parameters 
in our model, and so their combination lacks the flexibility 
to fit perfectly all data within the overlap region. This situa- 
tion could be improved by using additional free parameters, 
for example parameters controlling the rates of exponential 
decline, or by use of free-form fits (c.f. DP90). 

To compare the numbers of high-redshift objects in the 
samples with the models, Fig. ^ plots the un-normalized cu- 
mulative redshift distributions at z > 1.8 on a logarithmic 
scale. For the 3CRR sample, model B (that with no de- 
cline in p at high-redshifts) predicts four sources at z > 2 
(Qm = 1; for Qm = the model predicts three sources) 
whereas there is only one observed. Using Poisson statistics 
the probability of observing one or less objects given a mean 
of four is 0.1. Therefore the observations are marginally in- 
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Figure 9. Cumulative redshift distributions at z > 1.8. Each 
panel shows the distributions predicted by each of models A, 
B & C as in Fig. The data for each sample is shown as a 
solid line with the filled region the ±lcr Poisson uncertainty. The 
top three panels refer to the 3CRR, 6CE and 7CRS samples, 
respectively, for S!m = 1. The next three panels show the corre- 
sponding plots for flyi = 0. The final two panels show the 7CRS 
cumulative redshift distributions including only sources with lu- 
minosities in the range 27.0 < log 10 L151 < 28.0 (Qm = 1) and 



27.4 < lot 



10 ^151 



< 28.4 



0). 



compatible with the model. However, we note that the lack 
of z > 2 objects in the 3CRR sample could also be due 
to a steepening of the RLF at log 10 L151 ~ 28.5, due to 
a maximum radio luminosity of AGN, perhaps related to 
the maximum black hole mass in galaxies at these redshifts. 
This could be included in RLF models, but we do not do 
this in this paper, since there is little constraining data at 
such high luminosities and we wish to keep the numbers of 
free parameters in our models to a minimum. Note that the 
new equatorial sample of powerful radio sources defined by 
Best, Rottgering & Lehnert (1999) which has selection cri- 
teria producing a sample very similar to the 3CRR sample 
does indeed contain four sources at 2 > 2, so the lack of 



such sources in the 3CRR sample may well be due to small 
number statistics. The evolution of the RLF for such high- 
luminosity sources is discussed in more detail in Jarvis et al. 
(2000). 

Moving on to consider the 6CE and 7CRS redshift dis- 
tributions (Fig. []) , we find that all three models are consis- 
tent with the high-z data to within la uncertainties. There 
is no evidence here for a decline in the co-moving space den- 
sity of radio sources at high-redshift. 

The main benefit of the addition of the 7CRS com- 
plete sample to the existing 3CRR and 6CE samples, is 
that it enables the evolution in the luminosity band 27.0 < 
log 10 L151 < 28.0 (Qm = 1) to be determined over the wide 
redshift range z « 0.3 to z « 3 (see Fig. [j]). DP90 found 
stronger evidence for a redshift cut-off in the steep-spectrum 
population by considering only the sources within a simi- 
lar luminosity band (accounting for the difference in selec- 
tion frequency). Therefore we now repeat our analysis using 
only the sources within this luminosity band (note that for 

= 0, the band chosen is that corresponding to equiv- 
alent source flux at z = 2, i.e. 27.4 < log 10 L151 < 28.4). 
The bottom two panels of Fig. ^| repeat the 7CRS cumu- 
lative redshift distributions, but now including only those 
sources in these luminosity bands. For Qm = 1, the no cut- 
off model B now appears inconsistent with the z > 3 data. 
The model predicts four sources at z > 3 whereas only one is 
observed. This gives the same probability as we found with 
the z > 2 3CRR sources, i.e. a probability of 0.1 that the 
model and data are consistent. For the case of Qm = 0, only 
two sources are predicted at z > 3 by the model, so this is 
certainly consistent with the observation of one source. 

Thus we find very marginal evidence (90% confidence 
level) that the numbers of objects in the 27.0 < log 10 L151 < 
28.0 luminosity band falls slightly below the no cut-off 
model prediction. There is however, one systematic effect 
which might reduce even this significance. The expected 
numbers of high-redshift sources were calculated from the 
models assuming a spectral index between 151 MHz and 
151 x (1 + z) MHz of cz = 0.8. However, if high-z sources 
were systematically significantly steeper than a = 0.8, then 
the expected number of high-z sources would be smaller. 
This is because a source of a given luminosity on the flux- 
limit with a = 0.8, would clearly fall below the flux-limit 
if a were greater than this (see figure 1 in Blundell et al. 
1999). 

Our studies of complete samples (Blundell et al. 1999) 
do not reveal any intrinsic correlation between redshift and 
a evaluated at 151 MHz (0151), however the correlation be- 
tween luminosity and CV151 will mean that any high-redshift 
sources (which must be high luminosity) will tend to have 
steeper spectral indices. The consequences of this effect 
are discussed in detail in Section 3 of Jarvis & Rawlings 
(2000) and we defer interested readers to that paper. Re- 
calculation assuming a = 1.0 at high-redshifts for model 
B gives an expectation of only two sources at z > 3, 
thereby removing any discrepancy between this model and 
the data. Finally, we note that the chosen luminosity band 
of 27.0 < log 10 L151 < 28.0 excludes the four highest red- 
shift objects in the combined sample (see Fig. |l|) which have 
luminosities just above this limit. Extending the luminosity 
range considered to 27.0 < log 10 L151 < 28.5 would give a 
similar result to those found for the entire samples. 
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Figure 10. Banded V/Vmax test for the 356 sources in the 3CRR/6CE/7CRS complete samples. The left panel is for SIm = 1 and the 
right panel for SIm = 0. The error bars are equivalent to la. There are no significant differences between the form of the V/V m ax tests 
for the two cosmologies. 



4 THE V/Vmax TEST 

In this section the V/Vmax test is applied to the complete 
sample data. This test calculates the ratio of the volume V 
in the Universe enclosed by each source of redshift z to the 
maximum volume Umax, corresponding to the maximum red- 
shift z m ax a source of this luminosity would have if its flux- 
density lies above one or more of the survey flux-limits. If 
there is no cosmic evolution then sources would have values 
of V/Vmax randomly distributed between and 1 (Rowan- 
Robinson 1968; Schmidt 1968). Hence by calculating the val- 
ues of V/V ma:x for all the sources in a sample and averaging, 
one can determine the mean evolution of the sample. A value 
of < V/Vmax > greater than 0.5 corresponds to positive red- 
shift evolution and less than 0.5 indicates negative evolution. 
The distribution of < V/Vmax > expected in the case of no 
evolution is, for a large sample size N, a Gaussian centred 
on 0.5 with a = (127V)- 1/2 . 

Avni & Bahcall (1980) devised a method for combin- 
ing samples with various flux-limits and sky areas for the 
application of a generalised V/Vmax test. This version com- 
bines several samples such that they can be treated as a 
single sample with a variable sky area depending upon the 
flux-limit. The new test variables are the volume enclosed 
by sources V c and the volume available in the new sample 
V a . Note also that the fact that the 6CE complete sample 
has an upper as well as a lower flux-limit means that the 
volume available for a source in the 6CE sample is given by 
V a = Vmax — Vmm- To calculate values of V, one must make 
some assumption about the shapes of the radio spectra to 
determine values of z ma x. For this test it is assumed that 
radio spectra are not significantly curved over the region of 
interest, and spectral indices at a rest-frame frequency of 
151 MHz are used. 

The known strong positive evolution from z — to 
z ~ 2 would completely mask any change in the evolution 



at higher redshifts in the standard V/Vmax test. Therefore 
a banded version of the test is used here (e.g. Osmer & 
Smith 1980). This modification calculates the mean V e /Va 
restricting the analysis to z > zo at several different values 
of zo- Hence all the evolution below zo is masked out of the 
analysis. The statistic is now 



V.-Vo 

Va-Vo' 



(15) 



where Vo is the volume enclosed at redshift zo- 

Fig. [To| shows this banded V/Vmax test for the complete 
samples. In the banded test, the adjacent points are not sta- 
tistically independent and the error bars are much larger 
than the typical dispersion of points. Note that at low- 
redshifts (z < 1.5) the strong positive evolution of the radio 
source population is clearly seen with the values of 14/14 
significantly above the 0.5 no-evolution line. At z = 1.5 the 
values approach the 0.5 line where they remain (within the 
errors) out to the redshift limit of this sample (z = 4.0). 
The positive bump at z ~ 2.8 is marginally significant, but 
at this point there are very few sources left in the sample, 
so no firm conclusions can be drawn from this. Thus there 
is no evidence from the V/Vmax test for a high-2 decline (or 
increase) in the comoving space density of radio sources. 

The V/Vm&x test can also be used to study the evolution 
of the low- and high-luminosity populations separately. The 
complete samples were separated into two groups, depend- 
ing upon their emission line luminosities as described in Sec. 
3.4. These groups are supposed to signify the division be- 
tween objects in the low- and high-luminosity populations 
in the RLF. The banded V/Vnax test was performed for each 
group and the results shown in Fig. Both populations un- 
dergo positive evolution at low-redshift (z < 0.5) with the 
evolution stronger (and at a higher statistical significance) 
for the high-luminosity sources. At z > 1 most sources in 
the complete samples belong to the high-luminosity popu- 
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Figure 13. The evolving radio luminosity function of our model 
C (solid lines) plotted along with the PLE model of DP90 (dashed 
lines). The curves correspond to z = 0,0.5, 1,2 (bottom to top). 
The DP90 model has been transformed to 151 MHz by assuming 
all sources have spectral indices of 0.8. 



Figure 11. Banded V/Vm&x test for the high- (top panels) and 
low- (bottom panels) luminosity populations separately. The left 
panels are Om = 1 and the right panels S!m = 0. 




Figure 12. Banded V/V m ax test for the complete samples sep- 
arated into sources in different luminosity ranges. The ranges in 
log 1( , L151 included in each panel are labelled in the bottom-left 
corners (Qm = 1 only is shown here). 



lation, so its V/Vmax values are essentially the same as those 
in Fig. ^ for the whole dataset. There are hints of contin- 
uing positive evolution in the regime 0.5 < z < 1 for the 
low-luminosity sources (this is more apparent for £Im = 1 
than Q,m = 0). This evolution of the low-luminosity pop- 
ulation is consistent with that derived in the model RLFs, 
where the positive evolution ceases at z ~ 0.7. 



In Fig. |l^ we show the results of separating the sample 
into luminosity bands and repeating the V/Vm^ test. The 
< z < 1 evolution appears to depend strongly on luminos- 
ity. The luminosity separation does not provide any evidence 
for negative evolution at high-redshifts (i.e. a redshift cut- 
off) in either of the two high-luminosity bands. Therefore 
we conclude from this Section that the V/Vma,* test shows 
positive evolution at z < 2 which depends strongly on lumi- 
nosity, but is inconclusive at higher redshifts. 



5 COMPARISON WITH PREVIOUS WORK 
5.1 Comparison with DP90 

In Fig. [l3|, our model-C RLF is plotted with the steep- 
spectrum PLE model of DP90 at various redshifts (Qm = 1 
only). The DP90 LDE model which incorporates negative lu- 
minosity evolution at high-redshift is not shown here, since 
it is almost identical to the PLE model in the redshift range 
0.5 < z < 3 (which is to be expected as these regimes are 
constrained by data). To account for the fact that DP90 
evaluated the RLF at high-frequency (2.7 GHz), their lumi- 
nosity scale has been transformed to 151 MHz assuming a 
global spectral index of 0.8. The kinks in the DP90 model 
at log 10 (Li5i/ W Hz - sr~ ) ~ 28.2 are due to the use 
of a seven term power-law expansion for the non-evolving 
low-luminosity RLF component: there are no such powerful 
sources at low-redshift in their samples (or indeed in the 
3CRR sample) and hence no constraint from the data, so 
this feature is a numerical artefact. As described in Sec. 3.5 
our modelling procedure is also imperfect: the bumps in 
our estimations of the RLF around the breaks in the high- 
and low-luminosity populations may be artefacts of the re- 
stricted number of free parameters in our models. 

In general, the DP90 PLE model and model C agree 
reasonably well. At low luminosities the PLE model has a 
steeper slope of 0.69 than the slope of model C (0.54). Again, 
this slope is only well constrained at low-redshifts and the 
cause of a differing slope is probably due to the use here 
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of a more recent determination of the local radio luminos- 
ity function. The break luminosity appears at exactly the 
same value for the two models. Note how our use of two 
populations with differential evolution mimics the evolution 
of the break luminosity in the PLE case. Since we know 
that radio sources are short-lived with respect to the Hub- 
ble Time, pure luminosity evolution has no physical mean- 
ing and there is clearly not one population of radio sources 
which fades over cosmic time. In contrast, the two popula- 
tion model has some physical meaning and mimics PLE by 
having only density evolution for both populations. 

The high-luminosity power-law slopes for the DP90 
PLE model and our model C are very similar. However, at 
2 = 2 there is a clear excess of high-luminosity (log 10 L151 > 
27.5) sources in the DP90 model. At log 10 L151 = 28 the 
DP90 RLF is a factor of 2.1 higher than that of our model 
C. At 2 — 3 and log 10 L151 = 28 this factor has increased 
to 2.4. The effect of this is shown by a comparison of the 
expected high-redshift sources in the 6CE sample for RLF 
models with no redshift cut-off in Dunlop (1998) and in this 
paper (Fig. ^|). The plot by Dunlop shows a predicted excess 
at 2 > 2 over the number observed and this excess continues 
beyond the last data point at z w 3.4. Dunlop takes this as 
further evidence for the redshift cut-off. However, in Fig. ^ 
(second panel from top) we see that the 6CE redshift distri- 
bution is consistent with our model B at all redshifts. How 
can these two models, both of which claim to be no cut- 
off models, have such different predictions for the expected 
numbers of high-redshift sources? The answer is that all no 
cut-off models simply take the RLF at its peak at 2 ~ 2 and 
'freeze' the shape and normalization up to high-redshifts (at 
least to z ~ 5). Therefore the DP90 models naturally pre- 
dict more of these luminous sources at high redshifts because 
they have more of them at 2 ~ 2. 

The next question is why do DP90 find a higher density 
of powerful sources at 2 < 2 < 3 than found by our study? 
The answer is probably due to the fact that a very high 
fraction of the most distant sources in DP90 did not have 
reliable redshifts: inspection of fig. 9 of DP90 shows that 
only a quarter of their 2 > 2 steep-spectrum PSR sources 
had spectroscopic redshifts with the remainder estimated 
from the K — 2 diagram. Dunlop (1998) reports that further 
spectroscopy has shown the redshifts in DP90 to be typically 
over-estimated (most likely due to the positive correlation 
between radio luminosity and if-band luminosity observed 
by Eales et al. 1997). If these sources are actually at system- 
atically lower redshifts then this explains why DP90 observe 
an excess of 2 > 2 sources. Extrapolating this high source 
density up to 2 = 5 then predicts many sources in no cut- 
off models. Note that because few of their redshift estimates 
were actually in the region of supposed decline (at 2 > 3) the 
redshift over-estimation has the opposite effect to the naive 
expectation that this new information strengthens claims 
for the redshift cut-off. Thus we have an explanation as to 
why DP90 find stronger evidence for a steep-spectrum red- 
shift cut-off than we find with the 3CRR, 6CE and 7CRS 
datasets (having virtually complete spectroscopic redshifts) . 

DP90 argued that their V/Vmax test on the steep- 
spectrum population showed clear evidence for a redshift 
cut-off beyond 2 = 2. However we do not find a similar 
behaviour in our low-frequency selected samples, with ap- 
parently no evolution between 2 = 2 and 2 = 4. Inspection 
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Figure 14. The evolving radio luminosity function of model C 
(solid lines) plotted along with the quasar RLF model C from 
W98 (dashed lines). The quasar RLF has been multiplied by a 
factor of 2.5 to account for the fraction of luminous quasars com- 
pared to radio galaxies in low-frequency selected samples above 
the critical luminosity log 10 Li5i ?s 26.5 (Willott et al. 2000). 
Also plotted as dotted lines is the high— luminosity part of the 
model-C RLF. The curves correspond to 2 0, 0.5, 1, 3, 2 (bot- 
tom to top). 



of figs. 12 and 13 of DP90 show that their evidence of nega- 
tive density evolution at 2 > 2 is highly tentative: their data 
are only about la away from the null-hypothesis of no evo- 
lution. The banded V/V m &K test suffers significantly from 
large statistical errors when only a few sources are left in 
the samples at the highest redshifts, and this tends to limit 
its usefulness in practice. Combined with the uncertainty in 
the DP90 redshift estimates, it is difficult to see clear evi- 
dence for a redshift cut-off in their V/Vmax tests. 



5.2 Comparison with the radio-loud quasar RLF 

In line with unified schemes (e.g. Antonucci 1993), the high- 
luminosity part of the RLF derived here consists of radio 
galaxies and quasars which are identical except for the an- 
gle our line-of-sight makes with the jet axis. In Willott et 
al. (2000) we showed that the fraction of quasars in the 
high luminosity population is 0.4. Therefore one would ex- 
pect a very similar RLF to that of the quasars in W98, 
with an increase in number density of ~ 2.5. Fig. |l4| plots 
the RLF derived for all radio sources in this paper and the 
quasar RLF from W98 multiplied by a factor of 2.5. At high- 
redshifts the two RLFs appear very similar (there are very 
few low-redshift/luminosity quasars in the complete sam- 
ples - see Willott et al. 2000 for a discussion of this). The 
high-luminosity power-law slope for model C is 2.3, com- 
pared with 1.9 for the quasar RLF. The peak for the high- 
luminosity Schechter function is at approximately the same 
luminosity (log 10 (.Li5i) ~ 26.6) as the break required by the 
source counts for the quasar RLF. The quasar RLF of W98 
had a flat slope below this break, whereas the models of this 
paper have very rapid declines with decreasing luminosity. If 
the models described here are correct, then the quasar RLF 
should also have this Schechter function form. A survey at 
fainter flux-limits would probe the break region at high- 
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redshift and hopefully constrain the number of quasars (or 
BLRGs) at lower radio luminosities. A sample of quasars 
with a flux-limit of Sisi > 0.1 Jy has already been defined 
by Riley et al. (1999), however incompleteness arising from 
the bright optical limit of this sample (as discussed by W98) 
will complicate attempts to make quantitative estimates of 
the shape of the RLF at and below the break luminosity. 



6 DISCUSSION 

Using low-frequency selected radio samples with virtually 
complete redshift information, we have investigated the form 
of the radio luminosity function (RLF) and its evolution. 
Our results are generally in good agreement with previous 
work by DP90 (Figjl|) and W98 (Fig. [lj), and most dis- 
crepancies can be explained by limitations of the modelling 
methods forced on this and previous studies by the sparse 
sampling of the L — z plane afforded by available redshift 
surveys of radio sources. For the reasons outlined in Sec. 
most notably our virtually complete spectroscopic redshifts, 
we believe our estimates of the steep-spectrum RLF to be 
the most accurate yet obtained. 

We find that a simple dual-population model for the 
RLF fits the data well, requiring differential density evolu- 
tion (with z) for the two populations. This behaviour mimics 
the effects of a pure luminosity evolution (PLE) model but 
with a more plausible physical basis. The low-luminosity 
population can be associated with radio galaxies with weak 
emission lines, and includes sources with both FRI and FRII 
radio structures. The comoving space density p of the low- 
luminosity population rises by about one dex between z ~ 
and z ~ 1 but cannot yet be meaningfully constrained at 
higher redshifts. The high-luminosity population can be as- 
sociated with radio galaxies and quasars with strong emis- 
sion lines, and consists almost exclusively of sources with 
FRII radio structure. The comoving space density p of this 
population rises by nearly three dex between z ~ and 
z ~ 2 but cannot yet be meaningfully constrained at higher 
redshifts. 

Our results on the RLF mirror the situation seen in X- 
ray and optically-selected samples of AGN. Low luminosity 
objects exhibit a gradual [oc (1 + jz)~ 3 ' 5 from z ~ to z ~ 1] 
rise in p with z which matches the rise seen in the rate of 
global star formation (e.g. Boyle & Terlevich 1998; Frances- 
chini et al. 1999) , and also the rise in the galaxy merger rate 
(e.g. Le Fevre et al. 2000) . The p of high-luminosity objects 
rises much more dramatically [oc (1 + z)~ 5 ' 5 from z ~ 
to z ~ 2]. The similarities between the cosmic evolution of 
radio sources and these different types of objects have also 
been pointed out by Wall (1998) and Dunlop (1998). 

We have re-investigated the question of whether there is 
direct evidence for a high-redshift decline in the comoving 
space density p of steep-spectrum radio sources. Applying 
the V/Vmax test to the 3C/6CE/7CRS dataset provides no 
evidence for any such decline, which is a different result from 
that obtained by DP90. A fundamental limitation is that 
both our complete samples and those previously studied by 
DP90 include only a few sources with z > 2.5 (e.g. just 12 
in 3C/6CE/7CRS), so small number statistics are a huge 
problem in assessing the evidence for or against a redshift 
cut-off. This is considered in detail in Jarvis et al. (2000). 
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Figure 15. The flux-limits on the radio luminosity - redshift 
plane for low— frequency selected samples. The shaded band shows 
the luminosity range near the break in the RLF which contains 
70 % of the luminosity density at high-redshifts. 



What would be the best way of clearing up this un- 
certainty regarding the redshift cut-off? Larger samples are 
clearly needed, but in the case of the bright samples like 
3CRR, there is insufficient observable Universe to improve 
much on the data already available. Significant progress on 
fainter radio samples is however achievable. The choice of 
the flux-density level at which to concentrate efforts for the 
next large redshift survey is clearly crucial. We believe that 
the vital factor is being able to efficiently probe objects at 
or just above the break in the high-luminosity RLF at high- 
redshift. The reason for this is that integration of p x L151 
over the RLFs plotted in Fig. [] shows that about 70% of the 
luminosity density of the radio population at z ~ 2 — 3 lies 
in the narrow luminosity range 26 < log 10 L151 < 27.5. A fu- 
ture large redshift survey should be targeted at probing this 
population, which is seen only to z ~ 1 in the 7CRS. Fig. 
|l5| shows flux-limits in the L151 — z plane for various sam- 
ples with this luminosity range shaded. The plot illustrates 
that S151 ~ 0.1 Jy is therefore the natural choice for a new 
large redshift survey. Fainter surveys would be less useful 
as a probe of the total luminosity density at high-redshift: 
inspection of Fig. [IB] shows that they would become com- 
pletely dominated by the low-luminosity population, 

To investigate the likely content of a S151 > 0.1 Jy 
redshift survey, we have used models A, B and C RLFs to 
predict redshift distributions for a hypothetical 1000-source 
survey (Fig. |l^) . Note that all models assume that the evo- 
lution of the low-luminosity population freezes at z ~ 1. 
The first thing to notice is that the low-luminosity popula- 
tion dominates the total number of predicted sources in the 
case of Qm = 1, whereas in the £Im = case the two pop- 
ulations become approximately equal. This same result can 
be seen in the source counts plots of Fig. ^| and arises be- 
cause the low-luminosity RLF is currently constrained only 
at low-redshift, and the ratio of high- to low-redshift avail- 
able volume is much larger in a Iow-Om Universe. Thus, 
although a large fraction of the sources will be at z < 1, 
we predict that in a Iow-Om Universe in which evolution 
of the low-luminosity population freezes at z ~ 1, of order 
50 % of the sources will be in the high-luminosity popula- 
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Figure 16. The predicted redshift distribution of a sample selected at a radio flux-limit of $151 > 0.1 Jy (five times lower than the 
7CRS) for !1m = 1 (left) and Qm = (right). The low and high luminosity populations are plotted separately, each for the three models 
A, B and C. The vertical axis is the number of sources per unit redshift. The plots are normalised such that the total number of sources 
in the sample equals 1000. 



tion and lie at high-redshift, irrespective of the strength of 
the redshift cut-off. This should make determination of the 
evolution in the luminosity density a rather straightforward 
process. There should be no longer be a problem of small 
number statistics, at least at z ~ 3. Adopting model C, 7% 
(5%) of a hypothetical 1000-source redshift survey would lie 
at z > 3 for J1m = 1 (or Qyi = 0), whereas the same model 
predicts 3% (2%) for a survey like the 7CRS. 

The pronounced high-redshift tail in the case of 'no- 
cut-off' model B in Fig. [L6| should be easily refutable with 
a 1000-source redshift survey at £151 > 0.1 Jy. For the 
symmetric-decline model A (and Qm = 1) one would expect 
about 35 z > 3 sources in such a sample, whereas for model 
B we would expect 130. Indeed, Dunlop (1998) has already 
claimed to have ruled out the no-cut-off case using indirect 
constraints on the redshift distribution of a LBDS (Leiden- 
Berkeley Deep Survey) sample flux- limited at Si. 4 ~ 2 mjy, 
or S151 ~ 10 mjy. We note from Fig. [IB] that such a sam- 
ple will be overwhelmingly dominated by the low-luminosity 
population, and hence is too deep to provide optimal con- 
straint on the evolution of the luminosity density of the RLF. 
Although, as we shall argue below, indirect evidence favours 
a gentle decline in p at high-redshift, Dunlop's argument 
in favour of a redshift cut-off is not water-tight. It is cur- 
rently predicated on an extrapolation between two samples, 
the PSR and the LBDS sample, neither of which is close 
to spectroscopic completeness, and which are separated by 
two orders of magnitude in flux-density; at high-redshift the 
PSR and LBDS sources lie on either side of the range in L151 
which is key to delineating the evolution of the luminosity 
density of the radio source population with redshift. 

An alternative approach to obtaining complete redshift 
information for large samples is to refine the sample selection 
criteria to exclude the large numbers of low-redshift sources. 
This is being attempted by several groups using spectral 



index and/or linear size criteria (e.g. Bremer et al. 1998; 
Blundell et al. 1998; Rawlings et al. 1998; De Breuck et al. 
2000). However, the numbers of high-z sources which are 
missed in these samples are uncertain, so they will always 
struggle to prove the existence of a redshift cut-off, although 
they can place a lower limit on the space density (Jarvis et 
al. 2000). 

Another alternative approach is to study the evolution 
of flat-spectrum radio sources. Jarvis & Rawlings (2000) 
have recently critically re-evaluated evidence in this area 
(Peacock 1985; DP90; Shaver et al. 1996; Shaver et al. 1999) 
and concluded that constant-p models for the most luminous 
flat-spectrum sources can be ruled out at the ~ 2a level. 
They favour a decline in p for these objects by a factor ~ 3 
between z ~ 2.5 and z ~ 5, i.e. behaviour intermediate be- 
tween models B and C considered in this paper. However, 
they also point out that it is not yet clear how this decline 
should be interpreted in the context of the steep-spectrum 
population. 



7 CONCLUSIONS 

• We have used three low-frequency selected surveys with 
essentially complete spectroscopic redshift information, giv- 
ing unprecedented coverage of the radio luminosity versus 
2 plane, to derive the most accurate measurement of the 
steep-spectrum RLF yet made. 

• We find that a simple dual-population model for the 
RLF fits the data well, requiring differential density evolu- 
tion (with z) for the two populations. 

• The low-luminosity population can be associated with 
radio galaxies with weak emission lines, and includes sources 
with both FRI and FRII radio structures; its comoving space 
density p rises by about one dex between z ~ and z ~ 1 but 
cannot yet be meaningfully constrained at higher redshifts. 
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The high-luminosity population can be associated with ra- 
dio galaxies and quasars with strong emission lines, and con- 
sists almost exclusively of sources with FRII radio structure; 
its p rises by nearly three dex between z ~ and z ~ 2. 
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